Thermal entanglement and efficiency of the quantum Otto cycle for the su(l,l) 

Tavis-Cummings system 
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The influence of the dynamical Stark shift on the thermal entanglement and the efflciency of 
the quantum Otto cycle is studied for the su(l,l) Tavis-Cummings system. It is shown that the 
degree of the thermal entanglement becomes larger as the dynamical Stark shift increases. In 
contrast, the efHciency of the Otto cycle is degraded with an increase of the values of dynamical 
Stark shift. Expressions for the efficiency coefficient are derived. Using those expressions we identify 
the maximal efficiency of the quantum Otto cycle from the experimentally measured values of the 
dynamical Stark shift 

PACS numbers: 



I. INTRODUCTION 



The Jaynes and Cummings (JC) [l| model is a 
paradigm of quantum optics describing a two-level sys- 
tem coupled to a single mode of the radiation field. Due 
to its relative simplicity and relevance to applications 
in diverse areas of modern optics such as cavity quan- 
tum electrodynamics 0, Q and quantum state engineer- 
ing with Josephson junction devices Q , JC model is still 
of current interest. The JC model and its generaliza- 
tions 0-0] is particularly useful for studies in the field 
of quantum information processing Q. A key issue in 
this context is the concept of the entanglement and the 
correlations between quantum states . In recent years 
attention was devoted to so-called quantum engines [1- 
[Tsj and the relation between the thermodynamics and the 
entanglement [H, HH. In particular, two aspects are ad- 
dressed: The relation between the degree of the entangle- 
ment and the efficiency of quantum engines [15| and the 
role of thermal fidelity as an indicator for a thermal phase 
transition [l^ [l^ ■ Generally, studies of thermodynam- 
ical effects in nanosized systems are expected to reveal 
the connection between quantum physics and thermody- 
namics contributing to the understanding of fundamental 
physical problems, such as the Maxwell's demon and the 
universality of the second law of thermodynamics 
Our focus here is on the so-called quantum heat engine 
that produces work based on quantum effects. For in- 
stance, for a quantum cavity one considers the quantum 
Otto cycle. This cycle contains two quantum adiabatic 
and two quantum isochoric parts and therefore preserves 
the volume of the cavity. As a theoretical model we uti- 
lize the su(l,l) Tavis-Cummings (TC) model Q includ- 
ing the dynamical Stark shift (DSS) effect @. This model 
describes two-photon transitions between the ground and 
the excited state via an intermediate state. The interme- 
diate state can be eliminated from the equations of mo- 
tion ,18, on the cost of introducing a dynamical Stark 
shift 0] . Recently, the influence of DSS on the system's 
concurrence and fidelity was studied [20j. The main goal 



here is to study of the relation between the efficiency of 
the quantum Otto engine and DSS. We consider a quan- 
tum thermodynamical cycle for experimentally accessi- 
ble situations. Specifically we consider a scenario where 
a system, consisting of atoms and a quantum cavity, is 
connected to two thermal baths one with a high and the 
other with a low temperature. As shown below, there is 
a direct relation between DSS and the efficiency of the 
quantum Otto engine; for a given DSS one can predict 
the efficiency of the quantum Otto engine and visa versa, 
i.e. from the evaluated efficiency one may infer the values 
of DSS. 



II. THERMAL CONCURRENCE 

The Hamiltonian of two TC atoms placed in an ideal 
cavity reads Q 
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{S+ + S+)a^ + (5^ + ^2") {a+y 



(1) 



where g is the coupling constant between the atoms and 
the radiation field of the cavity, a and a+ arc the photon 
annihilation and creation operators. ^ is the strength 
of the Stark shift that results in an intensity depen- 
dent transition frequency. The atomic two-level sys- 
tems are described by the spin operators = ^az, 
= i (cTj; ± icTy) , where (Jx^y^z are Pauli operators. In 
([1]) we neglect the kinetic energy ^S-r of the center of mass 



2M 

of the atoms (with momentum P and mass M). This is 

2^ y 0] , where d is the atomic 



valid as long as « d 



dipolc moment, V is the volume of the cavity, Eq is the 
electric constant and n is the number of photons in the 
cavity. Due to the presence of the Stark shift term the 
eigenfrequency of the system is not constant anymore. 
Usually similar problems arise in non-stationary dynam- 
ics of nonlinear systems, where the frequency shift has 
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a crucial consequences (2l|. Here in a stationary prob- 
lem we also expect to see consequences of the frequency 
shift. As the Stark shift term changes the system's energy 
spectrum it should also affect the thermodynamic char- 
acteristics of the quantum cycle. In the limit of strong 
5 ^ C or weak <? <C C coupling regimes between radiation 
field of cavity and atomic subsystem, in the interaction 
representation Hamiltonian (1) splits in two commuting 
parts [2^. However, for moderate ( ^ g we expect to 
see strong correlation effects. As usual, the radiation 
filed is assumed to be prepared in a coherent state with 
the distribution function W„ = ^71- e~" , where a is the 
mean photon number n. We use the following notation 
for the relevant Hilbert space vectors (<? and e stand for 
the ground and excited state) 

\eg,n + 2), \ge,n + 2), \ee,n), \gg,n + 4). (2) 

In this basis and for the resonant case, i.e. for ujq = 2aj, 
we construct density operator of the system. Since work- 
ing substance of the Otto engine is formed by two TC 
atoms we are interested in the density operator of the 
atomic subsystem. Note that the field is the subsystem 
with a large number of degrees of freedom and is prepared 
in a coherent state. Therefore, the state of the field is not 
influenced by the atom-field coupling interaction. In par- 
ticular the state of the field can never be identified exactly 
due to its probabilistic nature. Consequently, for quanti- 
fying the entanglement between the atoms wc should av- 
erage and trace out the field states which leads to a loss 
of coherence between the atomic and the field subsys- 
tems. Nevertheless, the atomic states are still correlated 
and we expect to obtain a nonzero concurrence. In the 
language of density matrix, the coherence of the atomic 
states means non-vanishing off-diagonal matrix elements 
of the reduced density matrix. In addition to the con- 
nection between the engine efficiency and the thermal 
concurrence, further important issues to be discussed are 
the consequences of the system's transition form pure to 
the mixed state on the engine efficiency. Using basis vec- 
tors (2) after tracing field states for the reduced density 

operator of the atomic subsystem pa ~ Tr j ( ) we 
deduce 
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(3) 



The partition function reads 



Z = Tr{e-^") = 2 (^1 + cosh (/? + 462)^ . 



(4) 



Here a = , 6 = ga . The eigenvalues of the Hamilto- 
nian ([T]) are 
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Va^ + 462, E2 ^0,E3^ 0, E4 = Va^+462. 



Using (3), (4) wc follow Ref. [l^ and evaluate the con- 
currence C according to 



C = max (0, v i?i 



R4) 



where the square roots correspond to the eigenvalues of 
the matrix 

R ~ {ay (g) (Jy)p* {cTy ® cry)p- 

As we seek to inspect the influence of the Stark shift 
we consider an asymptotic case that corresponds to large 
values of the Stark shift g << In this case eqs. (3) 
and (4) simplify and we obtain 
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Z = 2(1 + cosh /?a). 
Furthermore, we obtain for the eigenvalues 
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R 



3,4 



Rl,2 > R3A 
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and therefore 



C = max(0, =0. 



Obviously a large Stark shift leads to a disentanglement. 
The same result we obtain in the case of a small Stark 
shift. Indeed, considering g » ^ we infer from eq. (3) 
and (4) that 



P = 



I \ ai a2 a2 \ 

0.2 0,2 I ai 



(8) 
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In what follows we introduce the following notations 
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( - 1 + cosh 2f3b) , 02 



2Z 



sinh 26/3, 



and 



Z = 2(1 + cosh 26/3). 
The matrix R and its eigenvalues read 
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7, and hence 



C = max(0,-2Vi?3 



(10) 



This Result is quite simple and evident. The value of 
DSS is the key factor since for both limiting cases the 
system is disentangled. To consider arbitrary values of 
DSS we introduce the following notations 



26 2g ^a^ 
— = -r^ V = Pa = — , 
a ^ 1 
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1 X cosh [y vT^ 



1 X sinh [j/vT^ 



1 + x^ 



(11) 



(9) and rewrite the density matrix in the form 
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xbi 



h 2:61 + 1(61-62) 
l + x6i + f(6i+62 

(61 - 62) - (61 ^ 
(61 - 62) - (61 



xbi (61-62) (61-62) \ 
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-62) (l + a;6i) x6i 
+ 62) a:;6i 1 + a;6i J 
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2(1 + cosh y \/T 



Consequently, the matrix R = (ctj, (E) ay^p*{cry (E) <Jy)p reads 



(l+2bix)(a;+2(bi-fc2)) 
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(13) 
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The eigenvalues of the Matrix (fT2|) and consequently 
the thermal concurrence is a function of two parameters: 
The first parameter is the relation between the cavity- 
atom coupling constant and the Stark shift x = The 

second one describes the temperature dependence of the 

f- 2 

concurrence y = Fig.l is deduced from the calcu- 

lations using (|13|) . From this figure we conclude that 
the dependence on the first parameter is more relevant 
for the thermal concurrence. Therefore, the relation be- 



tween the cavity-atom coupling constant and the Stark 
shift is the key parameter. In the limiting cases » g, 
and ^ << g, as was shown before analytically, the ther- 
mal concurrence is zero. The concurrence increases lin- 
early with the parameter g/^ up to the maximal value 
Cmax{g/£,) ~ 0.5. This fact can be explained analyti- 
cally. Considering in eqs. (3) and (4) g^ as a small 
parameter and a/3 as a large parameter and retaining only 
the first order terms we find that R matrix reduces to 
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2gye 2.g/^ 
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Consequently the expression for the concurrence 
C^ig/O = max{0, g/£,{^/2 — 1)) is in hne with the nu- 
merical result. In the opposite case, i.e. for < g, and 
(^gj^^yi ^ 1 we infer for the concurrence a square root 
decay 



FIG. 2: 1)^ — ^ B system is connected with hot reservoir Th, 
structure of energy levels is conserved En while populations 
are changed from Pi2 (i = 1,2,3,4) to Pn (i = 1,2,3,4)- 
Work performed is zero. Isochoric part 

2) B — > C Adiabatic part, level populations is conserved 
Pn {i = 1, 2, 3, 4). While energy levels are changed from 
En (i = 1,2,3, 4)-to £;,2(j = 1,2,3,4) 

3) C — ?■ -D is the reverse process of A — > B 

A) D A adiabatic part. Populations are conserved Pi2(i = 
1,2,3,4) energy levels changed Ei2{i = 1, 2, 3, 4)-from to 
En{i = 1,2,3,4) 



C{glCj = TOax(0, 
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' 26/3 



1), (14) 



which is also in a good agreement with the numerical 
results. Namely in the regime of weak coupling between 
atom and cavity concurrence is linear function olg/^. In 
the regime of strong coupling we have square root law for 
decay of concurrence. 



III. EFFICIENCY OF QUANTUM OTTO 
ENGINE AND THE DYNAMICAL STARK SHIFT 



the thermal reservoir changes the population of the levels 
Pij {i = l,2,3,4.j = 1,2) (isochoric part), while dur- 
ing the adiabatic parts the performed work changes the 



structure of the energy terms Ei. 



l,2,3,4.j = l,2) 



(du e to the change of the atom-cavity coupling constant 
[23j). For Pij the first index i enumerates the level pop- 
ulation probabilities, while the second index denotes the 
probabilities before or after the heat transfer from the 
thermostat to the system. For Eij the first index de- 
scribes the energy levels while the second defines the en- 
ergy spectrum before or after performing work (i.e., the 
performed work changes the structure of the energy spec- 
trum as was stated above.) The quantum Otto cycle is 
shown schematically in Fig. 2. Following the standard 
notations, e.g. as in [isj . we write 



Having defined the thermal concurrences for our sys- 
tem we focus now on the relation between the values of 
the Stark shift, the concurrence and the efficiency of the 
quantum Otto engine. A quantum Otto engine with two 
two-level atoms as a working substance was proposed and 
discussed in details in [l^ . Here we recall very briefly the 
main ideas: Two TC atoms are placed in an ideal loss- 
less cavity which is connected to two thermal baths with 
different temperatures 13l = 13 h = The system 
performs a quantum Otto cycte consisting of two adia- 
batic and two isochoric parts. The heat exchange with 
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(15) 
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Pi2 =cxp(+/3LA2)/Zif; 

P^2 = { ^1-: (16) 



Zl ■ 



P41 =exp(-/3iA2)/ZL. 



Here 



Ai = + 4,9?, A2 = a\ e + ^9la^ a't 



and 



Zh = 2(1 + cosh /3hAi), 



=2(l + cosh/3HA2). 



Using (14), (15) and the standard definitions of the ther- 
modynamic quantities such as the energy of the system 
E, and the transferred heat Q as well as the work per- 
formed A we find 



E -^E^jPij 
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4 4 
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= ^£;.i(/'a -/',:2) = 



2Ai 



/sinh/3LA2 sinh/?/fAi 
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(19) 
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/sinh/3//Ai sinh/3LA2\ A2QH 



" V y Ai 

Finally for the efficiency coefficient we obtain 



^ Qh-Ql ^ ^ _ A2 ^ ^ _ y^^^^(2o) 

Ql Ai + (51/e)' 

31 = 9 + Sg,g2 = g- 

For the small values of the Stark shift we arrive at the 
conclusions 



Vi^ «9) 



Qh - Ql Sg f Sg 



Qh g \ g 



while for large shift we deduce 

Qh - Ql _ {6g)g 



vii » 9) 



'.H 



e 



vi^»g) <vi^«g)- 

Efficiency of Otto engine r](£^), (20) has a maximum for 
zero values of the Stark shift 



Vi^^0) = l-g/ig + 6g), 



when thermal concurrence is zero. Therefore we can ar- 
gue that quantum correlations in the atomic subsystem 
hinders realization of quantum Otto cycle with maximal 
efficiency. 



IV. CONCLUSIONS 

The purpose of the present paper has been to study 
the influence of DSS on the thermal entanglement and 
the efficiency of the quantum Otto cycle. For this pur- 
pose we considered the su(l,l) Tavis-Cummings system 
and showed that the degree of the thermal entanglement 
follows the values of the DSS. In particular, for vanish- 
ing DSS, i.e. ^ = the system is disentangled, and the 
degree of the thermal entanglement increases with DSS. 
The system becomes disentangled again only in the limit 
£, ^ g, i.e. for a weak atom-cavity coupling. The effi- 
ciency of the quantum Otto cycle is maximal for a small 
values of DSS 77(^ « g) > 'r]{£^ » g), i.e. when sys- 
tem is disentangled. On the other hand, rj decreases as 
77(^) w l/^^- Using the asymptotic expressions for the 
efficiency coefficient (20), one may identify the maximal 
efficiently of the quantum Otto cycle from the experi- 
mentally measured DSS. For ^ << g, the efficiency is 
77 ~ Sg/g, where Sg is the amendment of the atom cavity 
coupling constant due to the performed work. If ^ > > (7 
then we obtain for the efficiency 77 ^ {Sg)g/£,^- For a 
given value of the DSS, one can quantify the thermal 
concurrence as well. Again the key point is the relation 
^/g. For S,/g < 1, the concurrence is given in an explicit 
analytical expression C{g/£,) see Eq. (14). For the oppo- 
site case S,/g > 1 we refer to the expression above. In the 
strong coupling limit ^ << g or large DSS ^ >> g system 
is disentangled. From the physical point of view this re- 
sult is expectable. In both cases the Hamiltonian of the 
system (1) can be presented as the sum of two commut- 
ing parts H,nt = Hi{sl,sf) + i?2(si,sj), [i?i,^2] = 0. 
Therefore, the reduced density matrix of the atomic sub- 
system is a separable pa = pi ffi P2- In the language of 
thermal entanglement both limits are identical (entangle- 
ment is zero), while the engine works with the maximal 
efficiency only for ^ << Therefore, we can conclude 
that a zero entanglement is essential but not sufficient a 
criteria to reach the maximal efficiency. Comparing Eq. 
(3) with the Eq. (5) we see that with the increase of 
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DSS the system performs a transitfon from a pure state 
to a mixed state where the non diagonal elements of the 
density matrix are zero Eq.(5). The system produces 
maximal efficiency for ^ « g, being in the pure state 
Eq. (8). Therefore, we can conclude that for a maxi- 
mal efficiency the working substance should be in a pure 



coherent state and the transition to the mixed state de- 
grades the efficiency of the quantum engine. 
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